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ABSTRACT 

The central star of the planetary nebula (CSPN) M 2-29 shows an extraordinary R Coronae Borealis-like fading event in its optical 
lightcurve. The only other CSPN to show these events are CPD-56°8032 (Hen 3-1333) and V651 Mon (NGC 2346). Dust cloud 
formation in the line of sight appears responsible but the exact triggering mechanism is not well understood. Understanding how 
planetary nebulae (PNe) trigger dust obscuration events may help understand the same process in a wide range of objects including 
Population-I WC9 stars, symbiotic stars and perhaps Asymptotic Giant Branch (AGB) stars with long secondary periods (LSPs). 
A binary scenario involving an eccentric, wide companion that triggers dust formation via interaction at periastron is a potential 
explanation that has been suggested for LSP variables. Model fits to the lightcurves of CPD-56°8032 and M 2-29 show the dust 
forms in excess of 70 AU at the inner edge of a dust disk. In the case of CPD-56°8032 this radius is far too large to coincide 
with a binary companion trigger, although a binary may have been responsible for the formation of the dust disk. We find no direct 
evidence to support previous claims of binarity in M 2-29 either from the OGLE lightcurve or deep medium-resolution VLT FLAMES 
spectroscopy of the CSPN. We classify the CSPN as Of(H) with T sff = 50 ± 10 kK and log g = 4.0 ± 0.3. We find a mean distance 
of 7.4 ±1.8 kpc to M 2-29 at which the M v = -0.9 mag CSPN could potentially hide a subgiant luminosity or fainter companion. A 
companion would help explain the multiple similarities with D'-type symbiotic stars whose outer nebulae are thought to be bona-fide 
PNe. The 7.4 kpc distance, oxygen abundance of 8.3 dex and Galactic coordinates (£ = 4.0, b = -3.0) prove that M 2-29 is a Galactic 
Bulge PN and not a Halo PN as commonly misconceived. 
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1. Introduction 

Long-term photometric monitoring surveys such as MACHO 
(Alcock et al. 1997) and OGLE (Udalski 2009) have opened a 
remarkable window into the decadal variability of rare, short- 
lived phases of stellar evolution. Perhaps least understood are 
those variables that exhibit fading events due to dust obscura- 
tion. These include R Coronae Borealis stars (RCB, Clayton 
1996), Population-I WC9 stars (Veen et al. 1998; Kato et al. 
2002), symbiotic stars (Mikolajewska et al. 1999; Mennickent 
et al. 2008; Gromadzki et al. 2009), PNe (Mendez, Gathier & 
Niemela 1982; Cohen et al. 2002; Hajduk et al. 2008, hereafter 
HZG08) and perhaps AGB stars with long secondary periods or 

* Based on observations made with the Very Large Telescope at 
Paranal Observatory under program ID 079.D-07 64(A). 



LSPs (Wood et al. 1999, 2004; Soszyhski 2007; Nicholls et al. 
2009; Wood & Nicholls 2009). 

In all these systems but the LSP variables, it is reasonably 
well established that ongoing dust cloud formation in the line 
of sight is responsible for the fading events. The exact mecha- 
nism that triggers the dust formation is not always known but 
shock waves are believed to play an important role (Fleischer 
et al. 1992; Goeres & Sedlmayr 1992, hereafter GS92; Veen et 
al. 1998). Stellar pulsations or interactions with a binary com- 
panion can help produce these shock waves. There are no bi- 
nary RCB stars known which leaves pulsations as the only pos- 
sible factor (e.g. Lawson et al. 1990). Amongst Population-I WC 
stars the non-WC9 types form amorphous carbon dust via wind- 
wind collision (WWC) in binary systems (Tuthill, Monnier & 
Danchi 1999; Williams 2008). The WC9 stars appear to be dif- 
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ferent in that while the majority do form C-based dust, it's un- 
clear whether this is due to WWC or shocks associated with 
the clumpy winds of single stars. Ground-based observations of 
WC9 stars show that most do not vary above cr=5-10 mmag 
(Fahed, Moffat & Bonanos 2009) which strongly suggests pul- 
sations are not responsible for dust formation in these stars. 
Symbiotic stars are by definition binaries and pulsations are 
common in the evolved components (e.g. Mikolajewska 2001). 
Pulsations are known to occur mostly in H-deficient PNe (e.g. 
Ciardullo & Bond 1996; Gonzalez Perez, Solheim & Kamben 
2006), but our knowledge of binarity in PNe is still in its in- 
fancy especially for orbital periods ^ 1 day (Miszalski et al. 2009, 
2010b; De Marco 2009). 

The origin of LSPs in AGB stars has a long history of de- 
bate but no satisfactory explanation for them exists (Nicholls et 
al. 2009). Dust is clearly associated with the phenomenon but 
the presence of dust alone cannot distinguish between competing 
scenarios (Wood & Nicholls 2009). An alternative way to view 
the problem is to examine PNe as the progeny of AGB stars with 
LSPs and to look for evidence that would support a given sce- 
nario. Of all the candidate scenarios a wide binary model (e.g. 
Wood et al. 1999) would be the most straightforward to verify 
with the safe assumption that the companion survives the PN 
phase. In this scenario a very low mass companion in a large 
eccentric orbit interacts at periastron triggering dust formation. 
The dust formed may remain clumpy or it could settle into a sta- 
ble disk structure. If a disk is present prior to PN ejection the dust 
geometry would conspire to produce a bipolar or at least asym- 
metric nebula (Bond & Livio 1990), while the disk itself would 
not necessarily survive the ejection event. This relationship with 
asymmetric PNe (see also Wood et al. 2004), and particularly 
the strong resemblance between fading events in PNe and LSPs, 
makes PNe a potentially valuable tool to address the longstand- 
ing LSP problem. 

M 2-29 (PN G004.0-03.0) is the least studied amongst the 
three PNe known to show dust obscuration events, partly be- 
cause this behaviour was only discovered recently (HZG08). 
M 2-29 is infamous for its apparently low oxygen abundance of 
7.4±0. 1 dex which lies in the metal poor regime of Galactic Halo 
PNe (Pena, Torres-Peimbert & Ruiz 1991, 1992, hereafter P91, 
P92; Howard et al. 1997). Torres-Peimbert et al. (1997, hereafter 
TP97) found M 2-29 to have a high density (n e > 10 6 cm 4 ) 
circumstellar nebula in HST observations and consequently re- 
marked upon the spurious nature of ground based abundance 
measurements. We obtained medium resolution integral field 
spectroscopy of M 2-29 with the VLT to measure its chemical 
abundances prior to the discovery of the fading event. Miszalski 
et al. (2009) confirmed the dense circumstellar nebula of TP97 
and further suggested that it may indicate the nucleus is a sym- 
biotic star. Here we expand upon this claim within the context of 
dust obscuration events in PNe. 

This paper is structured as follows. Section [2] reviews dust 
obscuration in V651 Mon and CPD-56°8032. Section [3] dis- 
cusses the photometric variations of M 2-29 with a particular 
focus on whether there is any evidence for binarity. Section 
[4] presents the first deep, medium-resolution, spatially -resolved 
optical spectroscopy of M 2-29. Classification of the CSPN of 
M 2-29 is discussed in Sect.[5]and the properties of its nebula are 
covered in Sect. [6] The dust properties of M 2-29 are discussed 
in Sect. [7] Section[8]outlines the similarities between M 2-29 and 
symbiotic stars, and we conclude in Sect. [9] 



2. Dust obscuration in NGC 2346 and Hen 3-1333 

Apart from M 2-29 only V651 Mon of NGC 2346 (PN 
G215.6+03.6; Mendez et al. 1982) and the [WC10] nucleus 
CPD-56°8032 of Hen 3-1333 (PN G332.9-09.9; Pollacco et al. 
1992; Jones et al. 1999; Cohen et al. 2002) were known to show 
dust obscuration events. If there is a common binary explanation 
for the fading events in AGB stars with LSPs and these PNe, 
then it is fundamental to establish the binary status or otherwise 
of these PNe. A key question is whether clearly defined or iso- 
lated fading events repeat periodically presumably due to the in- 
fluence of a companion. An additional tool that may be used to 
infer binarity is the detection of a (circumbinary) dust disk with 
Keplerian kinematics (Van Winckel 2003; Bujarrabal et al. 2005; 
Deroo et al. 2007). Binarity appears to be a prerequisite for disk 
formation but much further observational evidence is required to 
test this hypothesis (Van Winckel et al. 2009; Chesneau 2010). 

V651 Mon consists of an A5V star and an invisible hot 
sub-dwarf in a 16 day period seen in both photometric and ra- 
dial velocity observations (Mendez & Niemela 1981; Mendez 
et al. 1982). Undoubtedly the longest monitored CSPN, V651 
Mon has been observed since the 1970s (see e.g. Schaefer 1983, 
1985; Kohoutek 1995) and is also monitored by the American 
Association of Variable Star Observers (AAVSCQ) and the All 
Sky Automated Survey (ASAS, Pojmanski 2002). Isolated min- 
ima as well studied as the one observed by Kato et al. (2001) 
are difficult to glean from the early literature, while the early 
AAVSO data show what appears to be a prolonged minimum up 
until HJD-2450000 = -3900 days. This could correspond to 
the quasi-stationary Phase C described by GS92 where an equi- 
librium persists between dust formation and radiative pressure. 
Table [Upresents the epochs and duration of recent isolated min- 
ima from Kato et al. (2001) and ASAS. It is inappropriate to 
claim any periodicity until further observations of minima are 
obtained. If a period of ~5 years is adopted, then a minimum 
should have been observed around HJD-2450000 ~ 1800 days, 
but it is possible the conditions for such an event were not suit- 
able for grain growth to occur. 

If fading events in V65 1 Mon were to be triggered by a wide 
binary it is clearly not either component of the observed 16 day 
binary. Unfortunately the radial velocities presented in the lit- 
erature are too few to confirm or deny such a companion. A 
circumbinary disk seems unlikely to surround V65 1 Mon since 
the A5V component of the binary CSPN is only slightly red- 
dened outside fading events (Mendez & Niemela 1981; Mendez 
et al. 1985). The evacuated inner zone or 'void' (Balick 2000) 
also means the survival of a disk after the PN phase is unlikely. 
Nevertheless, high-resolution (interferometric) observations of 
NGC 2346 are urged. 

CPD-56°8032 is the cool [WC10] nucleus of Hen 3-1333 
which De Marco & Crowther (1998) found to have T eff = 30 
kK, Voo = 225 km s and a high mass-loss rate of M = 4 x 10~ 6 
M Q yr _1 . A dust disk of inner radius 97+11 AU was resolved 
by the VLTI (De Marco et al. 2002b; Chesneau et al. 2006) 
and dual-dust chemistry is present (Cohen et al. 2002). FigureQ] 
shows the visual lightcurve of CPD-56°8032 observed by one of 
us (AFJ) that updates the previous version published by Cohen 
et al. (2002). The minima are tabulated in Tab. Q] and, as for 
V651 Mon, we suspect CPD-56°8032 is in a prolonged quasi- 
stationary state of minimum after HJD-2450000 = 3500 days. 

Cohen et al. (2002) described a number of plausible sce- 
narios to describe the minima which have largely remained 
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Table 1. Isolated minima in V651 Mon and CPD-56°8032. 
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untested. Phenomenologically the lightcurve of CPD-56°8032 
resembles the events in RCB stars but over much longer 
timescales (e.g. Alcock et al. 2001; Tisserand et al. 2008). To 
gain further insight into the problem we can model the lightcurve 
using either the RCB star model of GS92 or the WC9 model of 
Veen et al. (1998). The main free parameters in these models are 
the radius Rd between the star and the forming dust cloud, and 
the dust outflow velocity i.e. the terminal velocity of the stellar 
wind Voo • The WC9 model of Veen et al. (1998) did not match the 
curvature of the ingress for any Rd and was not considered fur- 
ther. For the ingress we used Eqn. 33a of GS92 which is a simple 
quadratic independent of both v OT and Rd that only describes dust 
grain growth. For the egress we used Eqn. 35 of GS92: 

R 2 

m(t) - m{t\ ) = (m(to) - m{t\ )) — d -—r ( 1 ) 

(R d + Voo t) 2 

where m(to) and m{t\) are the magnitudes at light minimum and 
maximum, respectively. A remarkably good fit was found using 
the empirically determined Rd - 91 + \ \ AU and = 225 km 
s- 1 (Fig. [Q). 

The quality of the fit suggests the same physics apply to fad- 
ing events in PNe as RCB stars despite the very different time 
scales involved. In other words, it is likely that the dust forms at 
the inner edge of the dust disk in CPD—56°8032. As Rj is 3—4 
orders of magnitude larger than in RCB stars, presumably be- 
cause of the stronger ionising field of the CSPN, the wind den- 
sities would be insufficient to support dust formation (De Marco 
et al. 2002a). The only way to reach the high densities required 
would be either with a dust disk or via wind-wind interaction 
with a companion. The dust disk is already established and the 
latter scenario is unlikely in CPD-56°8032 because for an as- 
sumed orbital period of ~5 years, a primary mass of 0.6 M© and 
a main-sequence companion of 0.4 M©, the orbital separation, 
and therefore Rd, would be at most 13 + 2 AU. This radius is far 
too small to explain the fading events in CPD-56°8032 . 

In summary, there is insufficient evidence for obscuration 
events in PNe to be triggered by a wide interacting companion. 
Firstly, there is a lack of a strong periodic repeatability from both 
V651 Mon and CPD-56°8032 , and secondly, we have excluded 
a companion to CPD-56°8032 in a P ~ 5 year orbit as the trig- 
ger is based on the substantially larger dust formation radius at 
97 ±11 AU. While binarity may be responsible for the formation 
of a dust disk, binarity itself does not seem to play a strong role 
in the dust formation events. All the above suggests the events 
are triggered by stochastic instabilities in the dust disk itself, per- 
haps as a result of clumpy, high mass-loss rate winds from the 
central stars. Further observations of V65 1 Mon to clearly estab- 
lish the presence of a dust disk are urged. 

3. The OGLE lightcurve of M 2-29 

Figure [2] depicts the OGLE /-band lightcurve after final pro- 
cessing to place all OGLE data on a uniform photometric scale 
(Udalski et al. 2008). The MACHO Cousins R lightcurve is also 




-2000 2000 4000 

HJD-2450000 (days) 

Fig.l. An updated version of the CPD-56°8032 lightcurve 
published by Cohen et al. (2002), black points, and ASAS 
(Pojmanski 2002), red points. GS92 model fits for R d = 97 ± 11 
AU and v x = 225 km s _1 (blue lines) to the lightcurve have 
(left to right) min/max magnitudes of 12.0/11.0, 12.2/11.1 and 
12.4/1 1.0, and last 1450, 1900 and 1 100 days. 
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Fig. 2. OGLE /-band (black points) and MACHO Cousins R 
(green points, shifted -0.9 mag) lightcurves of M 2-29. The 
blue lines are GS92 RCB star model fits. The inset shows the 
ingress in the old (circles) and final (points) OGLE-I lightcurves. 
Spectroscopic observations are arrowed and ranges analysed for 
periodic variability marked. The dashed vertical line demarcates 
OGLE-II from OGLE-III data. 



shown to indicate the duration of the deep minimum (Alcock 
et al. 1997; Lutz et al. 2010)0 In the following we examine 
the large-scale and short-scale /-band variations, as well as the 
OGLE V-band data, to search for evidence of binarity. 

3.1. Large-scale variations 

Zebrun (1998) first classified the CSPN of M 2-29 as an RCB 
star candidate, but it was not until HZG08 that the connection 
was made with the PN. As for CPD-56°8032, we found the fad- 
ing event of M 2-29 could be successfully modelled using the 



2 The correct MACHO FTS identification number is 101.21568.184 
and not 101.21568.639 as listed by Lutz et al. (2010). 
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simple RCB star model of GS92. Our model fit in Fig.|2]using 
Eqn.Q] assumes an upper limit of v ro = 1200 km s (Sect. [5TTb 
and adopts Rj - 70 ± 10 AU, / = 16.18 at light minimum and 
/ = 14.75 at maximum. The ingress model assumes total occul- 
tation of the star which may explain the difficulty in fitting the 
preceeding slow decline. 

The simple assumptions of the GS92 model are by Ockham's 
Razor preferable over more complicated alternative scenarios 
that involve occultation by a large dust disk (e.g. HZG08). The 
larger v m creates a more rapid recovery than in CPD-56°8032 
and the very large dust formation radius strongly suggests a dust 
disk is present in M 2-29 (see also Sect. |7J. The double-peaked 
[O III] profile detected by Gesicki et al. (2010) may be evidence 
for Keplerian rotation in the dust disk following the example of 
EGB 6 (Su et al. 2010). An alternative view may be the profile 
reflects a small bipolar outflow emerging from the disk. 

The slow ~0.4 mag variations at light maximum in M 2- 
29 are not seen in V651 Mon or CPD-56°8032. Only varia- 
tions reflecting the 16 day orbital period are seen in V651 Mon. 
The timescales are again much greater than in RCB stars where 
they are pulsation-driven. As we find no evidence for pulsations 
(Sect. 13.21 ), we speculate that they are the product of a contin- 
uously variable (perhaps clumpy) mass-loss rate. Further time- 
resolved spectroscopy of the CSPN to track changes in stellar 
parameters should clarify their origin. At this stage it is highly 
premature to assign any type of periodicity to the fading events 
given only one is observed. HZG08 claimed to interpret the 'dip' 
around HJD - 2450000 = 3200 days (i.e. the third OGLE-III 
season) to represent a 'secondary eclipse' from which they de- 
rived an orbital period of ~18 years. There are a number of prob- 
lems with this interpretation. Firstly, the nomenclature of stellar 
eclipses cannot be applied to fading events, and secondly, there 
is no evidence for any declining trend in the adjacent photome- 
try near the gap in coverage that would suggest a fading event 
had occurred. The latter point is supported by the photometry of 
the third OGLE-III season which monotonically increases rather 
than decreases towards earlier dates. 

3.2. Small-scale variations 

HZG08 also claimed to find a short period of ~23 days dur- 
ing the ingress seen in Fig. [2] Unfortunately this part of the 
M 2-29 lightcurve is very undersampled and HZG08 presented 
no quantitative evidence for periodicity. We detrended the same 
data and failed to find any significant periodicity in a Lomb pe- 
riodogram (Fig. 0). The strongest but statistically insignificant 
peak at 31.2 days in Fig. [3] is highly likely to be a sampling 
alias of order one month. A large component of the scatter ap- 
pears to originate from the pipeline as there are large changes 
between the original and reprocessed lightcurves (inset Fig. 
This further reduces the reliability of the data in this region. 
Low-level nebular contamination also needs to be considered 
for any contribution as artificial variability can be imprinted 
on non-stellar sources in the stellar-optimised OGLE reduction 
pipeline (Wozniak 2000; Miszalski et al. 2009). Miszalski et 
al. (2009) showed lightcurves of some particularly bright PNe 
which showed artificial periods of one year. In M 2-29 the I- 
band nebular contamination is weak enough (Fig. |4j to cause no 
problems for the pipeline but it is certainly not absent (HZG08). 

The purported appearance of periodic variability only dur- 
ing ingress by HZG08 is also inconsistent with the litghtcurve. 
If this were the case we would expect variations to arise only 
from pauses in dust formation as in some WC9 stars (Veen et al. 
1998). Otherwise, variability due to pulsations or a close com- 
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Fig. 3. Lomb periodogram of the detrended ingress into the fad- 
ing event. Horizontal lines are the formal significance levels and 
the arrowed peak at 31.2 days is a likely sampling alias. 




Fig. 4. OGLE-III /-band (left) and V-band (right) images of 
M 2-29. The same logarithmic scaling is applied in each 30x30 
arcsec 2 image where North is up and East is to the left. The con- 
toured /-band inset highlights the elongated nebula oriented at a 
position angle (PA) of 168° that was missed by HZG08. 



panion would typically be expected to also occur in the optically 
thin portion of the lightcurve at light maximum. If a V65 1 Mon- 
like binary were present in M 2-29 we should be able to see this 
variability. Two well-sampled regions at light maximum contain- 
ing 199 and 456 observations were examined for periodic vari- 
ability (Fig. Second order polynomials were fitted to each 
series and subtracted before rejoining the data for period anal- 
ysis. Figure [5] shows no significant periodic variability present 
at light maximum. We also rule out the presence of pulsations 
above <x = 0.02 mag with timescales of a few days or less which 
are somewhat common amongst Of(H) nuclei (Sec. 15. II Handler 
1998, 2003). Although the MACHO data are of lesser quality 
compared to OGLE, we nonetheless searched for periods up to 
1000 days at light minimum but found no significant variability 
(Fig.|5j. 

We conclude that M 2-29 does not have a close irradiated 
binary companion with P ^ 1 day (e.g. Miszalski et al. 2009), 
nor does it have a close V651 Mon-like companion with P ~ 2 
weeks (e.g. Mendez et al. 1982). 
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Fig. 5. Same as Fig.EJbut for the detrended OGLE-II and OGLE- 
III data at light maximum (left column) and the MACHO V 
(blue) and R (red) data during light minimum (right column). 
All marked peaks in the OGLE data are sampling aliases of 6 or 
12 months, and 1, 2 or 3 days. These peaks are absent from the 
less frequently sampled MACHO data. 



3.3. OGLE V-band data 

Having ruled out a close binary companion, there is still a possi- 
bility that a wider companion may be present. Sparser OGLE V- 
band data for M 2-29 are well-suited to look for colour changes 
due to a companion (Zebrun 1998; Pigulski et al. 1993; Fig. 2 
of HZG08). The severe nebula contamination (Fig. HJ leads us 
to suspect the pipeline processed data may be unreliable (Sects. 
3.3 and 3.4 of Miszalski et al. 2009). Using the spectroscopy 
of P91 and P92 we find V = 15.5 + 0.1 mag by comparing the 
flux at /15480 to the equivalent in Vega. Remarkably, this mag- 
nitude agrees well with the V-band magnitudes at light maxi- 
mum as remeasured by HZG08 despite the higher uncertainty 
with application of the difference image anaylsis method (Alard 
& Lupton 1998) to the heavily nebula contaminated V-band data 
(Fig.|4]i. This may be a result of the careful nebula subtraction by 
HZG08 and so we can trust their reprocessed V-band lightcurve. 
We must however assume that no deep minimum occured dur- 
ing the P91 and P92 epochs and this is reasonable given the large 
percentage of time spent at light maximum. 

The depths of the fading event are / ~1.3 mag and V ~ 2.1 
mag. A colour change also occurs in V - / from 0.7 + 0.1 mag 
before minimum to 1 .5 mag during minimum. Could this be due 
to a companion? Assuming this change is solely caused by the 
dust produced during the event gives A(V) = 2.1 mag. This cor- 
responds to A(I) ~ 1.3 assuming A^s ~ 8100 A inferred from 
the mean extinction curve of Cardelli et al. (1989). This value of 
A(I) is more suitable than A eS ~ 9000 A for the OGLE /-band 
since for V - / < 2 mag there is negligible difference between 
the OGLE and Landolt systems (Udalski et al. 2002). As A(7) 
matches the /-band depth of the fading event we cannot use the 
OGLE colours to infer the presence of a companion. Additional 
colour information or ideally NIR spectroscopy during the event 
would be required to investigate this further as the MACHO pho- 
tometry are of insufficient quality. 




-4 -2 2 4 

Aot (arcsec) 



Fig. 6. The HST F656N (Ha+[N II]) image of M 2-29 (HZG08) 
overlaid with approximate positions (dots) and footprints (lines) 
of the FLAMES mini-IFU for exposures in the grating settings 
LR01 (blue), LR02 (cyan), LR03 (green), and LR05, LR06 and 
LR07 (red). Differential atmospheric refraction shifts the 2 x 3" 
mini-IFU between each exposure and from the reported centre 
(white cross). The approximate mean seeing (FWHM) is shown 
by circles of diameters 0.68" (solid) and 0.84" (dashed) for non- 
LR03 and LR03 exposures, respectively. 



4. VLT FLAMES integral field unit spectroscopy of 
M2-29 

Spectroscopic observations of M 2-29 were made by BM 
and AA on 11 June 2007 under visitor mode program 
079.D-0764(A) with the FLAMES facility of the 8.2-m VLT 
UT2/Kueyen (Pasquini et al. 2002). Within the 25' diameter field 
of view FLAMES can deploy 15 miniature integral field units 
(mini-IFUs) to observe science targets, with another 15 reserved 
for sky measurements. Each mini-IFU consists of an array of 
20 square microlenses, each 0.52" a side, giving an effective 
rectangular aperture of ~2 x 3" with each corner being inactive. 
M 2-29 was observed with a mini-IFU centred at approximately 
a = 18 h 06 m 40.85 s , 5 = -26°54'55.8" that was placed 9.9' 
from the centre of the FLAMES pointing at a = 18 h 07 m 17.81 s , 
6 = -27°00'25.3". Coordinates for the centre of M 2-29 were 
chosen from the SuperCOSMOS Ho- Survey (SHS; Parker et al. 
2005) images which have a world coordinate system accurate 
to 0.3" rms. The non-configurable position angle (PA) of the 
mini-IFU was 165.84 degrees. At this orientation the mini-IFU 
favourably includes the central star and a substantial amount of 
the surrounding nebula emission (Fig.[6]l. 

Table [2] summarises the exposures taken at light maximum 
(Fig. |2]l with the GIRAFFE spectrograph equipped with the 
2kx4k EEV 'Bruce' CCD. We used six of the eight avail- 
able low-resolution grating settings: LR01, LR02, LR03, LR05, 
LR06 and LR07 to yield two continuous wavelength ranges of 
3697-5074 A and 5722-8314 A. The wavelength range, resolv- 
ing power R (A/AA) and exposure times made in each setting are 
given in Tab. [2] Longer exposure times in bluer settings targeted 
the faint temperature sensitive diagnostic lines [O III] /I4363 and 
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Table 2. Log of M 2-29 VLT FLAMES GIRAFFE observations. 



Grating setting 


A (A) 


R (A/AA) 


Exp. (s) 


UTC 


LR02 


3948-4567 


10200 


600 


01:1740 


LR02 


3948-4567 


10200 


3600 


01:2842 


LR05 


5722-6510 


11800 


1800 


02:29:44 


LR06 


6420-7168 


13700 


900 


03:01:03 


LR06 


6420-7168 


13700 


120 


03:17:00 


LR07 


7075-8314 


8900 


600 


03:19:57 


LR01 


3607-4085 


12800 


1800 


03:31:08 


LR03 


4486-5074 


12000 


1800 


04:02:19 


LR03 


4486-5074 


12000 


120 


04:33:17 



[N II] ,15755 required to measure high-quality chemical abun- 
dances (Sect. [63}. Shorter exposures were not considered as no 
lines were found to be saturated in the longer exposures. The 
simultaneous arc lamp was turned off except for spectrophoto- 
metric standard star exposures. The average seeing conditions 
were 0.68+0.06" (FWHM) excluding the 1800 s LR03 expo- 
sure (0.84+0.19"). The airmass ranged from 1.83 to 1.03 with 
changes during exposures limited to ~0. 1 or less, with the ex- 
ception of 0.39 for the 3600 s LR02 exposure. Sky transparency 
was generally poor with thin to thick clouds present (relative flux 
rms of 0.04-0.08). 

The science frames were reduced with the IRAF3 task 
dofibers after standard initial processing using ccdproc. We ef- 
fectively removed the electronic glow defect of the 'Bruce' CCD 
by including dark frames of the same duration in ccdproc All 
science frames were cleaned of cosmic ray events (CREs) with 
the excellent LA. Cosmic program (Van Dokkum et al. 2001). 
Aperture identifications were made in dofibers using fibre infor- 
mation provided in the OzPoz and FLAMES binary table ex- 
tensions with particular care taken to not misidentify or mis- 
place any apertures. Low order polynomial fits were made to 
an average fibre flat field that was later divided through all in- 
dividual flat field spectra. All science spectra were then divided 
by their corresponding normalised flat field spectrum. ThAr arc 
lamp spectra were also extracted and a wavelength solution was 
generated and applied to science spectra in the standard way us- 
ing arc lists sourced from the ESO FLAMES pipeline. 

Flux calibration was provided by reducing in a similar fash- 
ion the sub-dwarf O star Feige 66 (Oke 1990) taken with the 
same field plate earlier the same night (albeit with a different 
IFU). All 20 fibres were averaged to form a combined Feige 
66 spectrum from which the sensitivity function was derived. 
Spectra belonging to each grating setting were then flux cal- 
ibrated individually. Determination of a smooth spectrophoto- 
metric response across all individual settings is discussed further 
in Sect. [63] 

Figure [7] shows the spaxels selected for extraction of an av- 
erage CSPN (S) and outer nebula (N) spectrum. In each grating 
setting an average sky spectrum is subtracted which comes from 
15 dedicated sky fibres uniformly distributed in the 12.5' diam- 
eter FLAMES field-of-view. The selected CSPN spaxels were 
those four with the highest S/N in the continuum. There are mul- 
tiple components to the CSPN spectrum including photospheric 
absorption and emission lines, a dense circumstellar nebula and 
the surrounding larger nebula. Obtaining a representative nebula 
spectrum for chemical abundance measurements unaffected by 

3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. See http://iraf.noao.edu 



Table 3. M 2-29 CSPN properties. 



Property Value 

Spectral type Of(H) 

V 15.50 

T eB (kK) 50 ± 10 

log g (cm s~ 2 ) 4.0 ± 0.3 

M (M Q ) 0.60+jJ f 

Gravity distance (kpc) 6.9^q 

v 0O (kms- 1 ) < 1200 

M (M e yr 1 ) ~1Q- 6 



the dense stellar core is more difficult because of the small 2x3" 
mini-IFU footprint. An ideal position would be the SE corner of 
the mini-IFU, however differential atmospheric refraction pre- 
cludes coverage of blue wavelengths in this region. As a com- 
promise we have selected an average of two spaxels per grating 
setting to cover the region on or just beyond the northern rim. 

5. The central star of M 2-29 

5.1. Classification 

Figure [8] depicts our rectified spectrum of the CSPN with many 
common features identified. Overall the CSPN is a hot post-AGB 
star with He II absorption lines as first shown by P91. Our de- 
rived and estimated properties of the CSPN are summarised in 
Tab. [3] No features that could be assigned to a companion star are 
seen. A TMAF0NLTE model atmosphere analysis (Werner et al. 
2003; Rauch & Deetjen 2003) was applied to yield T eff = 50± 10 
kK from weak ionisation equilibrium of N IV and N V lines and 
a surface gravity of log §=4.0±0.3 cm s~ 2 from He II line wings. 
Our T e ff is consistent with the total absence of nebular He II 
emission and T e ff < 40 kK is unlikely as there are no He I ab- 
sorption lines. The model normalised mass fractions used were 
close to the solar values: H (7.413X10 1 ), He (2.505x10 '), C 
(2.169xl0~ 3 ), N (6.208xl0~ 4 ) and O (5.379xl0~ 3 ). An addi- 
tional H-deficient TMAP model was also attempted but no ad- 
equate fit could be reached. With log g known the distance to 
M 2-29 may be estimated following the procedure of Heber 
et al. (1984). We estimate the CSPN mass to be 0.60^-f m 
from the evolutionary tracks of Miller Bertolami & Althaus et 
al. (2006) and take the monochromatic flux at /15454 from our 
TMAP model. A distance of 6.9^^ kpc is found where the er- 
rors are dominated by the uncertainty in log g. 

Apart from the improved resolution of previously known 
stellar He II and C IV features (P91; P92), our deeper FLAMES 
spectra show many stellar emission lines of Si II, Si III, Si IV, 
N IV and N V that are seen for the first time (Fig. |5}. These 
features were identified as stellar because of their larger equiva- 
lent widths compared to the sharp, narrow nebular lines follow- 
ing Fig. 4 of Mendez et al. (1990). Winds are not incorporated 
into the model so not all emission lines are fit. The peculiarly 
narrow C IV emission lines for example are relatively common 
in hot, windy post-AGB stars including ~30% of the symbiotic 
stars studied by Mikolajewska, Acker & Stenholm (1997). The 
H-rich composition means that we can definitively rule out a 
[WR] classification. The narrow He II /14686 in emission with 
FWHM=3.7 ± 0.2 A and H-rich atmosphere classifies M 2-29 
as Of(H) under the well-defined scheme of Mendez et al. (1990) 
and Mendez (1991). M 2-29 meets all the requirements for an 
Of classification, i.e. He II /14686 FWHM < 4A, non-blueshifted 

4 http://astro.uni-tuebingen.de/~rauch/TMAP/TMAP.html 




Fig. 7. Selected spaxels averaged to form the CSPN (S) and outer nebula (N) spectra in each grating setting. Lighter colours indicate 
higher average intensity per spaxel reflecting the central star position. Contours are chosen to show the main nebula features (Fig. 
i). 



He II /14200, ,14540 in absorption and symmetric Hy absorption. 
An upper limit to Voo = 1200 km s _1 for M 2-29 may be inferred 
from the spectroscopically similar NGC 6826 (r e ff = 50+10 
kK, log g=3.9+0.2) which straddles the Of and higher Vco Of- 
WR types (Mendez et al. 1990; Pauldrach et al. 2004). 

5.2. The broad Ha profile 

Miszalski et al. (2009) first remarked upon the very broad Ha 
profile which has a full-width at zero intensity (FWZI) of -2000 
km s _1 (Fig.fTOb. Broad Ha profiles are relatively uncommon in 
PNe with only a few known including the aforementioned NGC 
6826 (e.g. Kudritzki et al. 1997, 2006; see also Kudritzki & Puis 
2000). Kudritzki et al. (1997, 2006) derived mass-loss rates of 
10~ 8 -10~ 7 M yr~' for these CSPN by fitting their Ha profiles. 
Larger mass-loss rates correspond to stronger Ha profiles. In the 
case of M 2-29 we therefore estimate the mass-loss rate to be 
~10~ 6 M yr _1 to produce the larger than normal Ha profile for 
PNe. Our estimate requires confirmation through detailed mod- 
elling. It is encouraging that Gesicki et al. (2010) found the same 
rate is required to reproduce the kinematic age of the inner neb- 
ula, but they did not identify the CSPN as the responsible mass- 
losing star. As there appears to be no luminous AGB star com- 
panion in M 2-29 the mass-losing star must be the Of(H) CSPN. 

Broad Ha profiles are more commonly observed in symbi- 
otic stars (Van Winckel et al. 1993; Munari & Zwitter 2002) 
where they originate near the hot primary (Quiroga et al. 2002; 
Skopal 2006 and ref. therein). Skopal (2006) studied the Ha pro- 
files of symbiotic stars during their active and quiescent phases 
to find mass-loss rates of 10~ 7 -10~ 6 M G yr _I and~10~ 8 M yr _1 , 
respectively. The geometry of the mass-loss was also modelled 
by Skopal (2006) and found to likely be in the form of an opti- 
cally thin bipolar outflow. Gesicki et al. (2010) suggested such 




Velocity (km/s) 



Fig. 10. Ha profile for CSPN (top) and outer nebula (bottom). 
The top profile has been offset vertically to match the lower one. 



a geometry for M 2-29 consistent with this result. Section|8]will 
discuss further similarities with symbiotic stars. 

5.3. Silicon emission lines 

There are three ionisation stages of Si visible in our deep 
FLAMES spectrum, e.g. Si II ,15041, ,15056, Si III ,14552, ,15739 
(Fig. |9]l and Si IV ,14654 (Fig. |9]). They are typically found in 
late [WC] type nuclei where they can reflect an extreme over- 
abundance of Si (Leuenhagen & Hamann 1998; Hultzsch et al. 
2007). Hultzsch et al. (2007) modelled the spectrum of M 1-37 
to find a Si overabundance of 28 times solar! There appears to 
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Fig. 9. FLAMES spectra of key regions with stellar emission lines and a diffuse interstellar band (DIB) labelled. Emission lines 
from the inner nebula of N III, C III, O II (see Fig.[8]and Mendez et al. 1990 for identifications) and [N II] /15755 are also present. 



be no explanation for this as Si is not expected to be enriched by disk might help explain the presence of the low ionisation Si II 

AGB nucleosynthesis (Werner & Herwig 2006). Leuenhagen & species around the 50 kK CSPN. 
Hamann (1998) suggested the overabundance may be evidence 
for additional mixing but could there be another explanation? 

The presence of a silicate-rich dust disk in M 2-29 (Sect.0 could 5A - Interstellar absorption features 

potentially explain this perceived overabundance. We speculate Mogt ab tion lines not fit b our XMAP model are diffuse 

that the St lines are formed via some interaction between the interstellar bands (Herbi 1995) The Na l D1 and D2 a bsorp- 

CSPN wind and the dust disk. Shielding supplied by the dust don ^ m resolyed Mq ^ onents each . Gesicki et aL 
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(2010) claimed a component at Vh e i = -100 km s _1 was caused 
by a superposition of the high velocity cloud (HVC) 002.9-06.2 
that lies 3.4° away from M 2-29 with V lsr = -88 km s _1 (Putman 
et al. 2002). The authors missed a closer HVC 004.8-06.2 (3.3° 
away) with Vi sr = +206 km s . As these clouds measure only 
~0.2x0.2 degrees across on the sky there can be no line-of-sight 
superposition. 

We measured equivalent widths (EW) of 0.82, 0.89 and 1.41 
A for the Na I D 1 lines, all of which lie in the saturated regime of 
the EW-E(B - V) relation Munari & Zwitter (1997). No estimate 
of E(B-V) can therefore be made from the Na I Dl lines. Munari 
& Zwitter (1997) offer K I ,17699 as an alternative indicator at 
higher values of E(B - V), however our low resolving power and 
shallow exposure time in LR07 prevent its use here. Distance 
estimates from the Na I lines are also impractical towards the 
Galactic Bulge given the complex nature of randomly distributed 
dust clouds (Koppen & Vergely 1998). 

6. Nebula properties of M 2-29 

6.1. Reddening and integrated fluxes 

The logarithmic extinction at HB, c, varies somewhat in the small 
amount of literature on M 2-29. Ground-based longslit spec- 
troscopy values include 0.82 (Acker et al. 1992), 0.97 (P91; P92) 
and 1.18 (HZG08). A value of 1.05 was measured from difficult 
to calibrate fibre spectroscopy (Exter et al. 2004), while the mul- 
tiple grating settings of our FLAMES spectra make them unsuit- 
able for measuring c. An independent measure of c is required to 
check these values given the strong stellar Ha emission from the 
CSPN that may have contaminated previous measurements. One 
method involves the use of suitable emission line fluxes of the 
whole nebula. This method includes an [O III] to HB conversion 
which requires the [O III] and Ha emission to be reasonably 
co-located. This is confirmed after inspection of the HST im- 
ages and helped by the small magnitude and variation of [N II] 
(,16548+,16583)/Ha across the nebula (0.05-0.2). 

Table E] lists integrated Ha, HB and [O III] ,15007 fluxes of 
M 2-29. Fluxes measured by Acker et al. (1992, A92) are longslit 
values scaled to the whole nebula. Other fluxes were measured 
directly from an [O III] imaging survey of Galactic Bulge PNe 
(Kovacevic et al. 201 1, Kl 1), the Southern Ha Sky Survey Atlas 
(SHASSA, Gaustad et al. 2001), and HST F502N and F656N 
images (HZG08). The HST fluxes as listed exclude the CSPN, 
while values in parentheses include the CSPN. The HB fluxes 
were calculated using [O III] A5QQ1/HB = 4.86 (an average of 
our FLAMES spectra) in all cases except Acker et al. (1992). 
The procedure for obtaining Ha fluxes from aperture photometry 
of SHASSA data is described by Frew, Parker & Russeil (2006). 

For the HST fluxes we obtained the science-ready pipeline- 
reduced images from the MAST archival. Aperture photome- 
try was performed with radii of 60 pixels (2.73") and 8 pix- 
els (0.36") chosen to cover the whole nebula (including CSPN 
and inner nebula) and the CSPN and inner nebula, respectively. 
Surrounding annuli of width 10 and 2 pixels were used to sub- 
tract the sky and nebula contributions in each case. Flux cal- 
culations followed the WFPC2 Photometry Cookbook that im- 
proves upon Holtzman et al. (1995). Known HST sensitivities 
and rectangular filter widths of 35.82 A (F502N) and 28.33 A 
(F656N) were used to find the final sky-subtracted flux in each 
aperture. Both SHASSA and HST Ha fluxes were adjusted by 

3 http://archive.stsci.edu/hst 

6 http://www.stsci.edu/hst/wfpc2/wfpc2_faq/wfpc2_nrw_phot_faq.html 



Table 4. Integrated fluxes of M 2-29 uncorrected for reddening. 



Line 




A92 


Kll 


SHASSA 


HST 


-log 




12.18 


12.14 




12.34 (12.27) 


-log 


F([0 III]) 


11.47 


11.45 




11.65 (11.59) 


-log 


F(Ha) 


11.46 




11.50 


11.58 (11.50) 



0.01 dex to remove the small (< 2%) contribution from [N II] 
emission in the filter bandpasses. The HST fluxes appear to be 
systematically fainter than ground-based measurements which 
may be explained by the higher spatial resolution. 

Since the same aperture was used in each image we consider 
the Ha to [O III] ratio to be reliable. We therefore find c = 0.95 + 
0.05 from the average HST derived value ofHa/HB = 5.82 (Tab. 
5). This value corresponds to E(B - V) — 0.65 or Ay - 2.0 
and is free of any contribution from the CSPN or inner nebula. 
Table |4] also explains the higher c = 1.18 measured by HZG08. 
Comparing the CSPN-inclusive Ha flux to the CSPN-free HB 
flux gives HafHB = 6.92 or c = 1.18. The stellar Ha emission 
is clearly distorting the value obtained by HZG08. The higher 
quality of the P91 spectra could explain the better agreement 
with their value of c = 0.97. 



6.2. Distance and ionised mass 

With our reddening of c = 0.95 + 0.05, measured angular diame- 
ter of 5.0 x 4.8 arcsec at 10% of peak intensity and an integrated 
Ha flux of -log F(Ha) = 11.58 (Tab. [4} we estimated the dis- 
tance to M 2-29 using the Ha surface brightness-radius relation 
(SBR) of Frew & Parker (2006). The average relation suitable 
for M 2-29 gives a distance of 7.7 + 2.3 kpc. We adopt a final 
distance of 7.4 +1.8 kpc which is a weighted mean of the grav- 
ity distance (Sect. lBTTT i and the SBR distance. M 2-29 is therefore 
a Galactic Bulge PN entirely consistent with its Galactic coordi- 
nates (PN G004.0-03.0) and small radius. 

The ionised mass of the nebula can now be calculated 
which can generally serve as a useful discriminant between PNe 
and symbiotic nebulae. Santander-Garcfa, Corradi & Mampaso 
(2007) note typical ionised masses of 10 4 -10 2 M Q for symbi- 
otic nebulae and 0.1-1.0 M for PNe. Frew & Parker (2010) 
give a range of 0.005-3 M Q derived from a larger sample of 
PNe. Symbiotic nebulae typically have smaller masses since the 
nebula forms from the tenuous wind of the evolved companion 
(see Corradi 2003 for further discussion), but this method is not 
foolproof because some symbiotic nebulae can reach 0. 1 M Q 
(Santander-Garcfa et al. 2008). Using the formula of Pottasch 
(1996): 

M neb (M Q ) = 4.03 x 10- 4 <E U2 d 5/2 HB l/2 8 V2 (2) 

where d is the distance (7.4+1 .8 kpc), Hp is the dereddened flux 
(0.46 x 10~ n erg ctrT 2 s _1 ), is the nebular radius (4.9"), and 
e = 0.3 is the volume filling factor. This gives M ne b — 0.08+0.05 
M Q which lies in between a symbiotic nebula and PN but is more 
consistent with a PN. The nebula does not resemble those around 
D-type symbiotic stars (e.g. Corradi 2003), however the multiple 
nebulae present are shared with D'-type symbiotics in which the 
outer nebula is presumably ejected by the white dwarf (see Sect. 
[8}. There is therefore no reason to suspect the nebula around 
M 2-29 is not a planetary nebula. 
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Table 5. Emission line intensities of the outer nebula 
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Fig. 12. Residual H I emission line intensities for Case B condi- 
tions and c = 0.95. The slightly larger scatter at A < 4000A is 
consistent with the lower intensities of these lines and the diffi- 
culty of instrumental and flux calibration at these wavelengths. 



6.3. Chemical abundances 

M 2-29 has long been considered a halo PN because of its low 
oxygen abundance of 7.4 + 0.1 dex derived from ground-based 
longslit spectroscopy (P91; Exter et al. 2004). This anoma- 
lously low abundance is entirely the result of low spatial reso- 
lution ground-based observations averaging over the strong spa- 
tial variation of [O III] ^4363/^5007 (TP97; Miszalski et al. 
2009). TP97 measured an O abundance of 8.5 dex from HST 
spectroscopy of the inner nebula and yet it is still incorrectly re- 
ferred to as a halo PN in the literature (Otsuka et al. 2009; Lutz 
et al. 2010; Stasihska et al. 2010; Schonberner et al. 2010). To 
settle the issue we have recalculated the chemical abundances of 
M 2-29 from our spectrum of the outer nebula. 

Figure QT] displays the joined FLAMES spectra of the outer 
nebula and the inner nebula. Note the dramatic difference be- 
tween the outer nebula and the core which shows much stronger 
[Ne III], He I and [O III] 4363 emission lines in addition to lines 
not seen in the outer nebula (e.g. [K IV] ,16102, [CI IV] ,17530, 
8045). Here we focus only on the outer nebula for which it is 
more straight-forward to calculate abundances. Even after flux 
calibration of the six individual spectra there are still some un- 
determined offsets between them. These offsets were found by 
fixing c = 0.95 (Sect. |67T1 > and fitting the H I and He I emission 
lines to their expected Case B intensities assuming a Howarth 
(1983) reddening law. Figure [T2l shows the results for H I and 
demonstrates the corrected spectrophotometric solution is suit- 
able for abundance determination. Table lists our corrected 
emission lines intensities before and after dereddening. 

The measured line intensities were analysed via the plasma 
diagnostic program HOPPLA (Acker et al. 1991; Koppen et al. 
1991; see also Girard et al. 2007). The extinction constant was 
fixed to c = 0.95 and an iterative process was used to con- 
sistently derive the electron temperatures of T e ([N II])=8980 
K, T e ([0 III])=10500 K and an electron density n e = 2220 
cm -3 . Application of the usual ionisation correction factors al- 
lowed for the derivation of the elemental abundances in the for- 
mat 12 + log(X/H) in Tab. [6] Table [6] includes for comparison 
abundance measurements of M 2-29 from P91, TP97, Exter et 
al. (2004, E04), as well as average abundances for Bulge PNe 
(Cuisinier et al. 2000) and Solar abundances (Asplund et al. 
2005). Overall the abundances are in good agreement with TP97 
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287.28 


C II 6578.0 


0.81 


0.40 


[N II] 6583.3 


84.53 


41.70 


He I 6678.2 


9.09 


4.35 


rc TTi en 1 c a 

[i> 11J 6716.4 


4.48 


2.12 


[S II] 6731.6 


7.46 


3.52 


He I 7065.3 


10.53 


4.51 


[Arlll] 7135.6 


27.62 


11.61 


He I 7281.4 


1.77 


0.72 


[O II] 7320.0 


4.62 


1.85 


[O II] 7330.0 


3.97 


1.59 


[Ar III] 7750.9 


3.95 


1.43 



except for neon. Our measured oxygen abundance of 8.3 dex is 
low but entirely consistent at the 2<x level with an average Bulge 
PN (Cuisinier et al. 2000). There is now no reason to consider 
M 2-29 as a Halo PN since the oxygen abundance is only 0.4 dex 
lower than Solar. 

7. Dust properties of M 2-29 

The large 70 AU radius at which dust forms around M 2-29 sug- 
gests the presence of a dust disk similar to CPD-56°8032 (Sect. 
[3} and this is also supported by NIR and MIR observations. 
Gesicki et al. (2010) showed an archival Spitzer IRS spectrum 
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Fig. 11. VLT FLAMES spectra of the outer (top) and inner (bottom) nebula. 
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Table 6. Chemical abundances of M 2-29. 



Element 


P91 


TP97 


E04 


this work 


Bulge 


Solar 


He 


10.97 


- 


11.00 


11.02 


10.90 


10.93 


N 


7.2 


7.9 


7.2 


7.8 


8.2 


7.8 


O 


7.3 


8.5 


7.5 


8.3 


8.7 


8.7 


Nc 


6.7 


8.0 


6.8 


7.4 




7.8 


S 


5.9 


6.5 


5.6 


6.6 


6.9 


7.1 


Ar 


5.3 




5.8 


6.0 


6.2 


6.2 


log(N/0) 


-0.1 


-0.6 


-0.3 


-0.5 


-0.5 


-0.9 


log(Ne/0) 


-0.6 


-0.5 


-0.7 


-0.9 




-0.9 


log(S/0) 


-1.4 


-2.0 


-1.9 


-1.7 


-1.8 


-1.6 


log(Ar/0) 


-2.0 




-1.7 


-2.3 


-2.5 


-2.5 




20 22 24- 26 2B 3D 32 34 36 

Wavelength (fan) 



Fig. 13. Archival Spitzer IRS spectrum of M 2-29. Crystalline 
silicate features are identified (see text) and the IRS spectrum at 
bluer wavelengths can be found in Gesicki et al. (2010). 

of M 2-29 covering 5.7-14.2 yum. A peak at 11.3 pm was identi- 
fied by the authors as a polycyclic aromatic hydrocarbon (PAH) 
band, but the absence of other bands at 6.2, 7.7 and 8.7 pm is in- 
consistent with this interpretation. The same feature is detected 
in the D'-type symbiotic HD 330036 by Roche, Allen & Aitken 
(1983) which they attributed to SiC dust. Figure [131 shows the 
IRS spectrum at longer wavelengths after basic pipeline calibra- 
tions and reductions were applied. To identify crystalline sili- 
cate features we compared the IRS spectrum of M 2-29 with the 
features identified in the ISO spectrum of NGC 6302 (Tab. 1 
of Molster et al. 2001) and in a set of evolved stars (Tab. 1 of 
Molster et al. 2002). We note the clear presence of the follow- 
ing minerals (features with S/N above 30): Forsterite (23.7, 27.6, 
31.2 and 33.6 pm), Diopside (25.0, 29.6 and 32.2 pm), Enstatite 
(24.5 and 34.9 pm) and a feature at 34. 1 pm described as a com- 
bination of Forsterite and Diopside. The silicate rich composi- 
tion alone is not enough for a dual-dust chemistry (Cohen et al. 
2002) which cannot be assigned to M 2-29 as there are no strong 
PAH bands at bluer MIR wavelengths. 

Table |7]contains available near-, mid- and far-infrared fluxes 
from many ground and space-based missions (Beichman et al. 
1988; Price et al. 2001; Skrutskie et al. 2006; Churchwell et 
al. 2009; Carey et al. 2009; Ishihara et al. 2010). The core of 
M 2-29 as detected by 2MASS has colours J - H = 0.41 and 
H - K s = 0.17, well within the region occupied by PNe and 
substantially below the area that characterizes the D-type sym- 
biotics (e.g. Corradi et al. 2008). The steepness of the optical- 
NIR spectral energy distribution excludes a major contribution 
by free-free emission from the wind of the hot CSPN. We have 
therefore subtracted the stellar photospheric contribution which 
we normalized in the optical. The Rayleigh-Jeans tail of this 
starlight effectively removes any possibility of a credible excess 
of NIR emission (i.e. from a stellar companion), even depending 
on 2MASS data alone. But we emphasise this does not rule out 
the existence of a companion of sub-giant luminosity or fainter. 



Table 7. NIR, MIR and FIR fluxes of M 2-29 uncorrected for 
reddening. 



Band 


A cm (pm) 


Flux (mJy) 


2MASS J 


1.23 


4.4 


2MASS H 


1.66 


4.1 


2MASS K s 


2.16 


7.9 


IRAC [3.6] 


3.56 


28.9 


IRAC [4.5] 


4.52 


40.5 


IRAC [5.8] 


5.73 


36.3 


IRAC [8.0] 


7.91 


71.2 


MSX A 


8.28 


151 


MSXC 


12.13 


227 


MSXD 


14.65 


216 


MSXE 


21.34 


186 


AKAR1 S18 


18.0 


1180 


MIPS [24] 


23.7 


1150 


IRAS [12] 


12.0 


850: 


IRAS [25] 


25.0 


2000 


IRAS [60] 


60.0 


2280 


IRAS [100] 


100.0 


<75000 



Longward of the K s -band there is a strong excess which we 
attribute to thermal emission by dust grains. We followed the 
dust fitting procedures common when handling dusty WCL Pop- 

I WR stars in the use of Colangeli et al. (1995) laboratory mea- 
surements of the wavelength dependence of amorphous carbon 
dust (e.g. Williams et al. 2001). Our results are shown in Fig. 
[14]for coherent combinations of the available photometry as de- 
scribed below. 

In the MIR there are two fortuitous combinations of photom- 
etry from two independent missions. MSX observed M 2-29 in 
1996-7 during the obscuration event, while Spitzer IRAC reob- 
served it in 2005-6 at light maximum. M 2-29 is not listed in the 
point source catalogues of MSX (PSC2.3) nor in the GLIMPSE- 

II source lists of Spitzer. Therefore, we measured M 2-29 in the 
absolutely calibrated MSX images in all four long wavelength 
bands (8.3-21.3 fim). Likewise we obtained photometry in the 
four IRAC bands (3.6-8.0 pm) to which we applied the appro- 
priate aperture corrections, and the MIPS 24 pm band whose 
diffuse absolute calibration is equivalent to that of MSX (Cohen 
2009). Fitting blackbody distributions to the data we find in the 
obscuration event 300K grains match all the MSX points, while 
outside of this both warm (550K) and cool (100K) grains are 
required to match the double peaked SED. 

The earliest relevant FIR data are those by IRAS taken during 
1983. Two dust components can be seen but from very different 
temperatures than those observed by Spitzer. IRAS requires 200 
K and 60 K grains to fit the 12-60 pm excess emission and, of 
course, it cannot observe the warm components seen by MSX 
and Spitzer IRAC. The IRAS 25 pm detection of 2.0 Jy (based 
on observations over a 6-month span) makes an interesting com- 
parison with both Spitzer MIPS 24 pm in 2005-6 (1.15 Jy) and 
AKARI's All-Sky Survey 19 pm flux of 1.18 Jy during 2006-7. 
All three FIR measurements were made at light maximum. The 
recent AKARI and Spitzer data are in excellent agreement and 
confirm a real change since IRAS, whose bright 25 pm observa- 
tion has an uncertainty of +1 1%. Note also that the temperature 
of the cool dust also cooled during this period. These changes are 
clearly consistent with the high mass-loss rate of an evaporating 
dust disk. 
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Fig. 14. MIR and FIR energy distributions measured by IRAS in 1983 (left panel), MSX in 1996-7 during the main obscuration event 
(middle panel) and Spitzer in 2005-6 (right panel). Blackbody fits of 200K/60K (IRAS), 300K (MSX) and 550K/1 10K (Spitzer) are 
also shown. 



8. Similarities with yellow symbiotic stars 

As there is currently no evidence for an evolved companion to 
the Of(H) CSPN the chance that M 2-29 is a symbiotic star 
is remote. However, it is still instructive to draw this compar- 
ison to demonstrate that there is indirect evidence for binarity 
in M 2-29. The closest sub-group of symbiotic stars to M 2-29 
are the yellow or D'-type symbiotics (Schmid & Nussbaumer 
1993). Examples include AS 201 (Kohoutek 1987), V471 Per, 
StHa 190 and HD 330036 (Schmid & Nussbaumer 1993). 
Instead of the usual M-type giant or Mira companion of D-types, 
D'-types have yellow G-type giants (Belczyhski et al. 2000). 
They also exhibit much cooler dust emission (200-400 K) than 
D-types (800-1000 K), a dense nebula of N e ~ 10 6 -10 7 ctrT 3 , 
and a relatively cool WD of T e $ ~ 50-60 kK. In this work we 
have demonstrated M 2-29 shares all these properties except for 
the yellow giant companion. 

The evolutionary status of D' symbiotics is discussed by 
Jorissen et al. (2005). Their WD components are thought to have 
recently evolved from the AGB which explains the remnant dust 
emission, relatively cool WD temperatures and ionised nebulae 
(Schwarz 1991; Van Winckel et al. 1994). AS 201 in particular 
matches the nebula configuration of M 2-29 with a dense inner 
nebula surrounded by an older planetary nebula thought to be 
ejected by the WD. Interaction with the wind of the cool giant 
companion creates the dense inner nebula. There is a small pos- 
sibility this is happening in M 2-29 if an obscured yellow sub- 
giant is present. A sub-giant would not yet be luminous enough 
to rival the luminosity of the Of(H) CSPN (M v = -0.9 mag 
at 7.4 kpc), especially if there is some additional obscuration 
from the dust. This would explain the its absence in our VLT 
FLAMES spectroscopy. Dusty cocoons appear in the symbiotic 
stars RX Pup (Mikolajewska et al. 1999), other symbiotic Miras 
(Gromadzki et al. 2009), the S-type yellow symbiotic AE Cir 
(Mennickent et al. 2008), and even in some single C-rich Miras 
like R Vol (Whitelock et al. 1997). All of these share with M 2- 
29 similar depths and durations of fading events due to dust. RX 
Pup is a particular interesting example where during obscuration 
all photospheric features of the Mira disappear (TiO bands in the 
optical and CO bands in the NIR), although pulsations can still 
be seen in the NIR. 

There are other explanations for dense inner nebulae in PNe 
but they all involve binary CSPN. Rodriguez et al. (2001) found 
a dense circumstellar nebula in Hen 2-428, a bipolar PN which 



has since been found to have a double-degenerate close binary 
CSPN (Santander-Garcfa et al. 2010). We discount this possi- 
bility for M 2-29 as we found no evidence for a close binary 
(Sect. [32). A more feasible scenario may be a very wide main- 
sequence companion as in EGB 6-like CSPN (Frew & Parker 
2010; Miszalski et al. 2010a). The archetype of this class ex- 
hibits a dense circumstellar nebula around a resolved late-M 
main-sequence companion and cool dust at T~300 K despite the 
highly evolved nature of the CSPN (Liebert et al. 1989; Fulbright 
& Liebert 1993; Bond 2009; Su et al. 2010). If this dust did 
not survive from the AGB phase, then the only explanation for 
its presence would be ongoing interactions with the companion. 
This may involve a mismatched interacting wind (Kenyon et al. 
1993; Nussbaumer 2000; Miszalski et al. 2010a), however this is 
highly uncertain given the limited understanding of wide binary 
interactions in general. It is interesting to speculate that EGB 6- 
like CSPN could be the final evolutionary stage of M 2-29. 

One way to unambiguously identify symbiotic stars is the 
detection of the Raman scattered O VI lines at -16825 and /17082 
(Schmid 1989). Unfortunately, the low T eff = 50-60 kK of M 2- 
29 and D'-type symbiotics (Schmid & Nussbaumer 1993) is not 
hot enough to produce O VI which is only seen in the hottest 
symbiotics that have enough scattering particles available (see 
Mikolajewska et al. 1997). The Raman feature is therefore a suf- 
ficient but not necessary condition for a symbiotic star classifi- 
cation. Furthermore, even if M 2-29 had a primary hot enough 
to excite O VI, the curious appearance of the Raman lines in 
AE Cir only during obscuration (Mennickent et al. 2008) means 
we might not expect to see them at all during light maximum 
when our VLT FLAMES spectra were taken. Another interest- 
ing example is CI Cygni which never shows the Raman lines and 
only episodically shows [Fe VII] lines (Mikolajewska & Ivison 
2001). Long-term spectroscopic monitoring of M 2-29 may sim- 
ilarly reveal evidence of a companion. 



9. Conclusions 

We have studied the three PNe whose CSPN are known to show 
dust obscuration events: Hen 3-1333, NGC 2346 and M 2-29. 
Analysis of the lightcurves suggests a binary companion is not 
responsible for triggering the events but may be responsible for 
establishing the dust disks surrounding their CSPN. A partic- 
ular focus was the unusual Galactic Bulge PN M 2-29 using 
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the definitive OGLE lightcurve, deep medium-resolution VLT 
FLAMES mini-IFU spectroscopy and infra-red observations. A 
number of misconceptions surrounding M 2-29 were clarified 
during the course of the study. Our main conclusions are as fol- 
lows: 

- An updated CPD-56°8032 lightcurve was modelled for the 
first time using the simple RCB-star lightcurve model of 
GS92. A good fit was obtained for a dust outflow velocity 
of Voo = 225 km s" 1 (De Marco & Crowther 1998) and a 
dust formation radius of 97 ± 1 1 AU (Chesneau et al. 2006). 
The dust causing the obscuration events forms at the inner 
edge of the dust disk resolved by the VLTI. There is no pos- 
sibility that a binary companion could trigger these events 
via interaction at periastron at such a large radius to produce 
the pseudo-periodic events in the lightcurve every ~ 5 years. 
If PNe with dust obscuration events are the progeny of AGB 
stars with LPVs, then this argues against binarity as the cause 
forLSPs. 

- A similar application to the OGLE lightcurve of M 2-29 finds 
a plausible inner radius for a dust disk of 70 + 10 AU. Mid- 
infrared spectroscopy shows the disk to be rich in crystalline 
silicates with no PAHs present. Recent AKARI and Spitzer 
observations at 24 pm confirm the flux has dropped from 2.0 
Jy from IRAS to 1 . 1 8 Jy reflecting the evaporation of the dust 
disk by the high mass-loss wind of the CSPN ~ 1 0~ 6 M yr~ 1 . 
Stochastic variations in the wind and/or dust disk are likely 
to be responsible for the irregularity in obscuration events. 

- Analysis of the OGLE lightcurve finds no evidence for bi- 
narity as previously claimed by HZG08. A colour change 
during the obscuration event is only caused by the dust and 
not a companion. 

- A TMAP NLTE model was fitted to the VLT FLAMES 
spectrum of the CSPN to find r eff = 50 ± 10 kK and log 
g = 4.0 + 0.3. We classified the CSPN as Of(H) following 
Mendez et al. (1990) and Mendez (1991). Many weak emis- 
sion lines were found in the spectrum including Si II, Si III 
and Si IV. These may not represent an overabundance of sil- 
icon, but rather some kind of interaction between the Of(H) 
wind and the silicate-rich dust disk. A gravity distance of 
6.9 + 3.0 kpc was found after assuming M = 0.60 M from 
evolutionary tracks (Miller Bertolami & Althaus 2006). 

- HST images of M 2-29 were used to independently derive 
the reddening to be c = 0.95 ± 0.05 or E(B - V) = 0.65 
mag. This value is free from the strong stellar Ha contami- 
nation that effects previous ground based measurements and 
was used in calculating a statistical distance of 7.7 + 2.3 kpc 
using the Ha surface brightness-radius relation of Frew & 
Parker (2006). A weighted average distance of 7.4 +1.8 kpc 
implies an ionised nebula mass of 0.08+0.05 M Q that is more 
consistent with PNe than symbiotic nebulae. 

- Chemical abundances of M 2-29 were calculated for a spec- 
trum of the outer nebula uncontaminated by the dense cir- 
cumstellar nebula that contaminated most previous studies. 
The abundances are in agreement with TP97 and the oxy- 
gen abundance of 8.3 dex, together with the 7.4 kpc distance 
and its Galactic coordinates (PN G004.0-03.0), clearly dis- 
proves the common misconception that M 2-29 is a Halo PN. 
M 2-29 is a Galactic Bulge PN. 

- M 2-29 shares a number of similarities with D'-type sym- 
biotic stars including a dense circumstellar nebula, an outer 
planetary nebula, a cool T e ff = 50 kK CSPN and cool dust 
(100-300 K). We can rule out a yellow giant companion but 
a subgiant luminosity class or fainter cannot be ruled out as it 



may be obscured by dust or hidden in the glare of the Of(H) 
My = -0.9 CSPN. Further observations to search for such a 
companion are required for a definitive answer. 
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